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Abstract. We consider initial-boundary value problems for the derivative non- 
linear Schrodinger (DNLS) equation on the half-line x > 0. In a previous work, 
we showed that the solution q(x, t) can be expressed in terms of the solution of a 
Riemann-Hilbert problem with jump condition specified by the initial and bound- 
ary values of q(x, t). However, for a well-posed problem, only part of the boundary 
values can be prescribed; the remaining boundary data cannot be independently 
specified, but are determined by the so-called global relation. In general, an effec- 
tive solution of the problem therefore requires solving the global relation. Here, we 
present the solution of the global relation in terms of the solution of a system of 
nonlinear integral equations. This also provides a construction of the Dirichlet-to- 
Neumann map for the DNLS equation on the half-line. 
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1. Introduction 
The derivative nonlinear Schrodinger (DNLS) equation 

iqt + Qxx = i{\q\ 2 q) x , (l-l) 

where q(x, t) is a complex-valued function, arises in the study of wave propagation 
in optical fibers [1J and in plasma physics [B] (see [7] for further references). It is 
an integrable equation and the initial value problem on the line can be analyzed 
by means of the Inverse Scattering Transform (1ST) as demonstrated by Kaup and 
Newell [8]. In the last fifteen years, a generalization of the 1ST to initial-boundary 
value (IB V) problems developed by Fokas and his collaborators [9J [10l [12] , has made it 



possible to analyze equations such as (1.1) on domains involving a boundary. Several 



of the most well-known integrable PDEs (such as the KdV, modified KdV, nonlinear 
Schrodinger (NLS), sine-Gordon, and Ernst equations) have been investigated using 
this approach, see e.g. [H El US [13 [HI UHl [20] . 
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Recently, the Fokas method was implemented to the DNLS equation (1.1) posed 
on the half- line x > Provided that the solution exists, it was shown in [TB] 

that the solution q(x, t) can be recovered from the initial and boundary data via the 
solution of a 2 x 2-matrix Riemanm-Hilbert (RH) problem. The jump matrix for this 
RH problem is given explicitly in terms of four spectral functions a(k), b(k), A(k), 
and B{k), where k G C is the spectral parameter of the associated Lax pair. The 
functions a(k) and b(k) are defined in terms of the initial data qo(x) = q(x,0) via a 
system of linear Volterra integral equations. The functions A(k) and B(k) are defined 
in terms of the boundary data go(t) = q(0,t) and gi(t) = q x (0,t) also via a system 
of linear Volterra integral equations. However, for a well-posed problem, only one of 
the functions go and g\ (or a combination of these two functions) can be specified; 
the remaining boundary data cannot be independently specified, but are determined 
by the so-called global relation. Thus, before the functions A(k) and B{k) can be 
constructed from the above linear integral equations, the global relation must first be 
used to eliminate the unknown boundary data. 

The analysis of the global relation can take place in two different domains: in the 
physical domain or in the spectral domain. Although these two domains are related 
by a transform, each viewpoint has its own advantages. 

In the first part of this paper, we analyze the global relation in the spectral domain. 
We will show that A(k) and B(k) can be determined via the solution of a system of 
nonlinear integral equations formulated explicitly in terms of the initial data and the 
known boundary values. 

In the second part of the paper, we analyze the global relation in the physical 
domain. We use a Gelfand-Levitan-Marchenko (GLM) representation to derive an 
expression for the generalized Dirichlet-to-Neumann map (i.e. the map which de- 
termines the unknown boundary values from the known ones). Once the unknown 
boundary values have been determined, A(k) and B(k) can be constructed from the 
linear integral equations mentioned earlier. This representation has the numerical 
advantage that the system of integral equations is defined on a bounded domain. 

In the case of the NLS equation, a construction of the Dirichlet-to-Neumann map 
was presented in P]lll]. The Dirichlet-to-Neumann map for the sine-Gordon equation 
as well as the two versions of mKdV were analyzed in [llj. The analysis of the 
KdV equation presents some novel difficulties which were finally overcome in |21j . 
Although our approach is inspired by the preceeding references, the analysis here 
presents a number of novelties: (a) The solution of the global relation presented 
in section [3] takes place entirely in the spectral domain. This is in contrast to the 
approach of [11] , in which it is necessary to derive a GLM representation for an 
appropriate eigenfunction of the Lax pair before an analogous result can be obtained. 
The introduction of a GLM representation amounts to passing from the spectral to 
the physical plane. Here we remain in the spectral plane throughout the derivation, 
which simplifies the arguments both on a practical and on a conceptual level. The 
observation that it is possible to solve the global relation directly in the spectral 
plane was first made for the NLS equation in [14] . (b) In the special, but important, 
case of vanishing initial data go = 0, we will show that a further analysis of the 
global relation can be used to simplify the integral equations derived in the physical 
plane considerably. In fact, the number of equations reduces from five to two in this 



Physically, this type of initial-boundary value problem arises naturally. For example, assuming 
that we can create or measure the waves at some fixed point x = in space, we arrive at an initial- 
boundary value problem on the half-line with initial data given by the initial wave profile for x > 
and boundary data provided by the measurements at x = 0. 
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case. The case of vanishing initial data is of special interest in applications, where 
it can sometimes be assumed that the physical field under consideration is initially 
at rest before exteriorly created waves enter the domain. This is also the relevant 
initial condition in the application of the numerical methods based on the so-called 
"exact nonreflecting boundary conditions" [22J. A similar simplification for the NLS 



equation was presented in |14j . (c) The Lax pair for equation (1.1) contains terms 
of higher order in the spectral parameter k than what is the case for the examples 
above. However, by utilizing additional symmetries, calculations can still be kept at 
a reasonable length, (d) As noted in [16], the definition of eigenfunctions with the 



appropriate asymptotics of the Lax pair of ( 1.1 ) naturally leads to the introduction of 



a certain closed two-form A. In order to close the system of integral equations which 
characterizes the solution of the global relation, the system must be supplemented by 
an additional (somewhat complicated) equation for the second component of A. 

The paper is organized as follows: In section [2j we recall the Lax pair formulation 
and the global relation associated with (1.1). In sections [3] and [5j the global relation 



is analyzed in the spectral and physical domains, respectively. In section |4| we derive 
a GLM representation for an appropriate eigenfunction of the Lax pair. 

2. A Lax pair and the global relation 



Equation (1.1) admits the Lax pair [8] 

f ^ x + ik 2 [a 3 ,^>] = U X *, 
\ y t + 2ik 4 [a 3 ,V] = U 2 V, 

where a 3 = diag(l,— 1), k G C = C U {00} is a spectral parameter, ^(x,t,k) is a 
2 x 2-matrix valued eigenfunction, and the 2 x 2-matrix valued functions Q(x, t) and 
{Uj(x,t, k)}1 are defined by 

Q = ( ° j* ) , U x = kQ, U 2 = 2k 3 Q - ik 2 Q 2 a 3 - ikQ x a 3 + kQ 3 . 



Following |16| . we define a transformed eigenfunction fi(x,t, k) by 

<£(x, t, k) = e iJ &°) Aa3 fi(x, t, k)e if &% Aa3 , (2.2) 
where A = A\dx + is the closed real-valued one-form 

A(x,t) = ^\q\ 2 dx+(^\q\ 4 - l -(q x q-qq x )^dt. (2.3) 

We introduce Q(x,t) and Qi(x,t) by 

where 03 acts on a 2 x 2 matrix A by a 3 A = [03, A], i.e. e a3 A = e (73 Ae~ (7S . The 
function [i satisfies the Lax pair 



I fj, x + ik 2 [a 3 ,[i] = Vifi, 
I fit + 2ik 4 [a 3 ,fi] = V 2 H, 

where the 2 x 2-matrix valued functions {Vj(x,t, k)} 2 are given by 

V 1 = kQ-iA 1 a 3 , V 2 = 2k 3 Q-ik 2 Q 2 a 3 -ikQ 1 a 3 + kQ 3 - iA 2 a 3 . (2.5) 
Let T > be some given final time; we will assume that T < 00. Assume that q(x,t) 



is a solution of ( 1.1 ) with sufficient smoothness in the domain {0<x<oo,0<i<T} 
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and with sufficient decay as x — > oo. We introduce three solutions {fij(x,t,k)}f of 
(2.4) as the solutions of the following linear Volterra integral equations: 

Hj(x,t, k) = 1+ f X,t e i[kH-'-^kHt'-t)]a 3 ^ Vlli .)^ i if j k )dx' + {V 2N ){x',t\ k)dt') , 

(2.6) 

where (xi,ti) = (0,T), (^2^2) = (0,0), (2:3, £3) = (oo,£), and I denotes the 2x2 
identity matrix. The eigenfunctions {/ij}? satisfy the symmetries 



fij(x,t,k) = <riHj(x,t, k)ai, Hj(x,t,k) = (T3Hj(x,t,-k)(T3, j = 1,2,3. (2.7) 

The first of these symmetries implies that we can introduce complex- valued functions 
a(k), b(k), A(k), B(k), *i(t,fc), $ 2 (t,k) by 



s(k)= m Km s(k) _ / A(k) B(k)\ 

V ; \b(k) a{k) y K ' \B{k) A(k) J' V ' 

Zl JfO.O) " 



^i(£,fe)e J (°.°) *2 (*,*:) 

where s(fc) and <S(fc) are defined by 

*(fc) = Ms(0, 0, k), S(k) = [e 2tk4T& ^ 2 (0, T, k)]~ x . (2.9) 

Since a(k) = a(—k), the zeros of a(k) always come in pairs; if kj is a zero, then so is 
— kj. 

The function fi 3 (x,0,k), and hence also s(k), can be constructed from the initial 
data qo(x) via the linear Volterra integral equation 

H 3 (x,Q,k) = I+ f ' e ik ^ x '- x)&3 {V lf M 3 )(x' \0,k)dx' '. (2.10) 
J 00 

Similarly, ^2(0, £, k), and hence also S(k), can be constructed from the boundary 
values go(t) and g\{t) via the linear Volterra integral equation 

M 2 (0,£,&) = I+ I' e 2 ' tk ^'~ t)a3 (V 2 iJ2)(0,t',k)dt'. (2.11) 



Defining the sets D^, j = 1, . . . , 4, by (see figure [TJ 

= {(j - l)vr/2 < arg(fc 2 ) < jvr/2} , 

it follows from (2.10)-( 2.11[ ) that a(k) and 6(/c) are analytic and bounded in D\ UD2, 
while A(k) and B{k) are entire functions which are bounded in D\[JD 3 . The functions 
3>i(£, fc) and ^(i, k) are entire functions of k which are bounded for k G Z?2 U -D4. 

2.1. The global relation. The eigenfunctions fi 2 and ^3 are related by 

fi 3 {x, t, k) = fi 2 (x, t, k)e-^ k2x+2k4t ^ s s(k), (2.12) 

where the first and second columns are valid for k £ D 3 U D4 and k 6 Z?i U D2, 
respectively. The relation obtained by evaluating (2.12) at (0, £) is called the global 
relation. It imposes a relation between the Dirichlet and Neumann boundary values 
of q(x,t). The Dirichlet-to-Neumann map is determined by solving this relation for 
the unknown boundary values. Defining functions c(£, k) and d(t, k) by 

2i f (0 ' () A 

e J(o,°) 1 
c(t,k) = ~QjCJc\ — (^ 3 ( ' t > ^)) 12 ' d(t,k) = -—^(/j, 3 (0,t,k)) 22 , 
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>- Rek 



Figure 1. The sets {Dj}f that decompose the complex k-pl 



ane. 



we can write the (12) and (22) entries of the global relation as 

c(t, k) = *!(t, k) + e 2if ^ A ^$ 2 (t,k)e- 4lkH , k G D 1 U D 2 , 

a{k) 



d(t, k) = Mt, k) + e *«*M*^ e -«* 

a(k) 



keDiUD 2 



(2.13a) 
(2.13b) 



We find from the definition of /i3 that the functions c(t, k) and d(t, k) are analytic 
and bounded for k £ D\ U D 2 away from the possible zeros of a{k). 



3. The solution of the global relation: spectral domain 

In this section we analyze the global relation in the spectral plane. For concrete- 
ness, we consider the Dirichlet problem; the treatment of the Neumann problem is 



analogous. The main result is theorem 3.1, which expresses the spectral functions 
A(k) and B{k) in terms of the Dirichlet boundary data and the initial data via the 
solution of a system of nonlinear integral equations. 

For a subset E C C of the complex fc-plane, we let <9X denote the boundary of S, 
oriented so that £ lies to the left of <9£. We let £>t = Dj n {Im k > 0}, j = 1, . . . , 4. 
Moreover, we let 7 denote the contour dDf and let 7' denote a contour obtained 
from 7 by deforming it so that it does not surround the possible zeros of a(k), see 
figure [2j 

Theorem 3.1. Let T < 00. Let qo(x), x > 0, be a function of Schwartz class and 
let go(t), < t < T, be a sufficiently smooth function. Assume that qo and go are 
compatible at x = t = in the sense that qo(0) = go (0). Let a{k) and b(k) be defined 
for k € D\ U D 2 by equations (2.S)-[2A^. Suppose that a(k) has a finite (possibly 
empty) set of simple zeros. We denote these zeros by {kj,—kj}^ and assume that 



Then the spectral functions A(k) and B(k) defined by (2.8)-(2.9) are given by 



A{k) = $ 2 (T, k) B{k) = -$i(T, k)e 



4ik 4 T -2i J T A 2 (0,t)dt 



(3.1) 
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Figure 2. The contour 7 in the complex k-plane and the deformed contour 7'. The 
zeros of a(k) are symmetrically distributed with respect to the origin. 



where the complex-valued functions $i(t,k), & 2 (t,k) and the real-valued function 
A 2 (0, t) satisfy the following system of nonlinear integral equations: 

'*i(t,fc)\ _ A) N 



+ 



t / e **4(f-t) [(-i\ go \2 k 2 + iA2 )$! + k {\g \ 2 g + 2k 2 g + i 9l ) <& 2 ] (f , fc) 
[fe (|ffo| 2 5o + 2/c 2 5o - %) $1 + (iM 2 *: 2 + *A 2 )$ 2 ] (t', k) 



1 



A 2 (0, t) = -Mt)| 4 - 4Im(g (t) C ( 3 )(t)) 



(3.2a) 
(3.2b) 



2\g (t)\ 2 Im 



k [$ 2 (t,k) - $ 2 (t,ik)]dk 



where the functions gi(t) and c^(t) are given for < t < T by 



(t)= Z -\go(t)\ 2 go(t) + 4c^(t) + 



91 



4«?o(i) 



7T 



k [$ 2 (t,/c) - $ 2 (t,ik)] dk, 

igo(ty 



k 



dk 



(3.3) 
(3.4) 



+ e 



2i/* A 2 (0,t')dt' 



2i r r M) - 



. k 2 ^ 2 (t-k)e-^ kH dk + A V kf^-^t^e- 4 ^ 
7r J a(k) ^ J a(kj) 

kj&D+ 



Proof. Letting gi(t) := q x (0,t), equations (3.1) and (3.2a) follow from the definition 
of $1 and <3? 2 together with (2.4). In order to derive (3.3), we substitute the expansion 



(1) (2) (3) 

Ha H) fi) 



1 



k^oo, j = 1,2, 3, 



(n) (n) 

into (2.4) and solve for the s. By the second symmetry in (2.7), y - is a diagonal 
(off-diagonal) matrix for n even (odd). We find (see [TB] for further details of this 
type of argument) 

1 



(i) 



-a 3 Q, 2ia 3l if = Qpf + |<2 3 + ^Qi, 
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and 



„ \QQx + ~QV 3) $ = -flat + ~-Q*Q,. + t O;,t :! 7 OO rj .,, 



(2) 



1 o2. 

4-" 



4 l 



Thus, as fc — >• oo, 



1 /•(*,*) 



2fe fc 3 l 4 



2 



(x,t) 



Co' 



M 2 g +o 



bj,j(x,t,k)) 



22 



1 + 



fc 2 



w + O 



(3.5) 



where the closed one-form u is defined by 



Co' 



M dx - - (\q\ 6 + 3i\q\ 2 qq x + |fe| 2 - gfe) dt. 



The expansions in (3.5) are valid for k approaching oo within the regions of bound- 
edness of (fJ>j)i 2 and (jij) 22 . Hence, 



c(t, fe) = 
*i(t,fc) 



2fc fc 3 U 5 



toft) , c (3) (*) , 



2k 



A; 3 



where 

c< 3 >(t) 
It follows that 



9l Wo 
4 2 



(o,t) 



-)■ 00, keD 1 UD 2 , 
, fc — > 00, k £ D 2 U Da, 
k ->■ 00, k £ D 2 UD^, 
(o,t) 



(3.6) 



Co' 



150 1 50, 



(0,0) 



*?>(*) 



Co'. 



(0,0) 



<7i(*) 



ffo| 2 5o + 4c( 3 ) + 2^ ^ 2) , 



< t < T. 



(3.7) 



The function <J>2(i, &) — & 2 (t,ik) is an entire function of fc and satisfies 

1 



$ 2 (t,k)-$ 2 (t,ik) = ~$ { 2 2 \t) + 



k -»• 00, k £ D 2 U Da. 



Letting Cr, > 0, denote the large arc of radius i? which closes .D^ at infinity, 

C R = {Re ia G C I vr/4 < a < vr/2}, 

we find 

k [* 2 (i, fc) - $ 2 (i, ifc)] dfc = -^^ 2) (t). 



A; [$ 2 (*,*0 - <Z> 2 (t,ik)]dk 



lim 

R— >oo 



Cr 



Consequently, 



*?>(*) 



2 



k [® 2 (t,k) - <S> 2 (t,ik)\ dk. 



(3.8) 



A change of variables k — > ik in the part of t he i ntegration that runs along the 
imaginary axis shows that the contour dD^ in (3.8) can be replaced by —7. Thus, 
equations (3.7) and (3.8) yield (3.3). 

Substituting the expression (3.7) for g\ into the definition (2.3) of A evaluated at 
(0, t), we find 



A a (0,t) = -\go\ 



4Im(5oc (3) ) -2|(7o| 2 Re^ j - 



(2) 
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In view of (3.8), this yields (3.2b). 

It remains to derive the expression (3.4) for c^ 3 \t). Assume first that a(k) has no 
zeros. Let e > 0. The function 



k 



c(t, k) + 



igo(t) 



l 



+ 



l 



k - ee- 7 ™/ 4 k + 



ee 



is analytic in D\ U D 2 and of 0(1/ k 2 ) as k — > 00. Thus, the sectionally analytic 
function 



k 



c(t,k) + 



jgojt) 



1 

k+ee—*'/ 4 



-ik 



jgojt) 



1 



c(t, ik) + m 2 ' [ ik _ € ^i/i + ik+ee -xi/4, 



k £ Di U D 2 , 

k G D 3 U At, 



satisfies a scalar RH problem with jump across R U iR. Let T denote the oriented 
contour 

T = [-zoo, 0] U [0, -00] U [too, 0] U [0, 00]. 

We find 



k 



c(t,k) + 



1 



+ 



1 



1 

2wi 



k - ee -7 ™/ 4 k + ee" 7 ™/ 4 
3 

k(M) + &(M0 + ^E 



igojt) 

4 ^Li s l- ee'™/ 4 



s=0 



dl 

T^k' 



(3.9) 



Thus, 



c (3) (t) + 



so(t) 



f 



k 



kc(t, k) + i/cc(t, ik) 



igo(t) 



Letting e — > and using the global relation (2.13a) to eliminate c(t,k) and c(t,ik) 
from the integrand, we find 

m{t)~ 



1 

2wi 



A: 



k 



dk 



(3.10) 



2vri 



6(fc) 



e- 4ife4 *dA;. 



a(/c) ' a(ik) 

The integrand in the first term on the right-hand side is analytic in D 2 and of 0(1/ A; 3 ) 
as k — > 00 in D 2 . Thus, the contour of integration V = dD\ + dD 2 in this term can 
be replaced by dD\. Using the change of variables k — > —ik, we find that the second 
term on the right-hand side of (3.10) equals 



1(0,0) 



TTt 



k' 



m 

a(k) 



-Aik A t 



dk. 



(3.11) 



The glob al rela tion (2.13a) together with the asymptotics (3.6) imply that the inte- 
grand in (3.11) is of 0(l/k 3 ) as k — > 00, k 6 D 2 . Thus, the contour of integration in 
(3.11) can be replaced by dD\. We arrive at 

c^(t) =--.'. / //- <!>,(/. A-') + /<I>i (/.//,■) + •••"'/•••'. elk 



1 

2vri 



(3.12) 



iri 



k 2h[ ^® 2 {t,k)e- AikH 



dDi 



a{k) 



dk. 
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The symmetry properties (2.7) imply that both integrands in (3.12) are odd functions 
of k. Thus, we obtain 



.(3) 



(*) 



I 

+ 



k 



*i(t,fc) (*,«'*) + 



(o,t) 



TT 



k 



dk 



a{k) 



dk. 



This proves (3.4) in the case when a(k) ^ 0. 

If a{k) has a finite number of simple zeros {kj, —kj}^ C D\ U D2, we replace 
equation (|3.9|) with 



k 



c(t, k) + 



1 



1 



ee" 



-TTi/4 



+ 



1 



lc(t,l) +ilc(t,il) + 



igo(t) 



k + ee" 7 ™/ 4 
3 



N 



E c ^) 



3=1 



1 



+ 



1 



s=0 



ee 



-7Tl/4 



E c ;C)E 



where Cj(t), j = 1, 



s=0 

, iV, denotes the residue of c(t, fc) at fc = fcj, i.e. 



dl 
l-k' 

(3.13) 



Ci(t) 



By closing the contour at infinity, we find 

A 7 " 



1 

2^i 



3 

,s=0 



i s l 



dl 



l-k 



N 2 N ( 1 

2 E^ + pE c A 2 + ° U4 



as k — > 00 within Di U -D2- The contribution to from the term on the right- 
hand side of (3.13) that involves the residues {cj}^ is therefore 2^- =1 CjAy. The 
effect of the term involving the Cj's on the left-hand side of (3.13) is to add another 
2^- =1 CjA;| toc^ 3 ^(t). Finally, when replacing the contour in (3.11) withS-Di, we now 
pick up a sum of residues from the zeros of a(k) in D2; this contributes an additional 
— 4^fe. £ £)+ Cj/c| to c^ 3 ). The total net contribution to c^ 3 ^(t) from the residues at the 
zeros of a(k) is therefore 

N N 

2 E c ^ 2 + 2 E c ^ 2 - 4 E <**? = 4 E c A 2 > 

i =1 i =1 %ecj fej6-D+ 



which proves (3.4) also when a(k) 7^ 0. 



□ 



In [IB], the solution q(x,t) was presented in terms of the solution of a matrix RH 
problem with jump condition specified by the spectral functions a(k), b(k), A(k), and 
B(k). Substitution of the expressions (3.3) and (3.4) into (3.2) yields a system of 
nonlinear integral equations involving the functions $i(i, A;), $2 and A2(0, t). 
Assuming that this system has a unique solution, A(k) and B(k) can be determined 
from (3.1 ). 

When a(k) 7^ 0, the functions A{k) and B(k) only enter the formulation of the 
RH problem as the combination A(k)/B(k) for k E dD%. Then it is sufficient to 
consider the system satisfied by k), $2(t,k), and A2(0, t) for < t < T and 
k G dDiUdDs; the function A(k)/B(k), k 6 dD%, can be recovered from the solution 
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of this system. In the case when a(k) has zeros, the residue conditions also need to 
be taken into account. 



Remark 3.2. 1. By replacing in the right-hand side of ( |3.4[ ) the contour 7 by the 
contour 7' which does not surround the zeros of a(k) (see figure [2]), the sum over 
kj 6 Df can be absorbed into the preceding integral. 



2. The integrals along 7 which appear on the right-hand sides of ( 3.2b )-( 3.4) are 



in general, not absolutely convergent; they are defined as the limits as R —¥ 00 of the 
integrals obtained by replacing 7 with 7^, where 7^ = [-Re 7 ™/ 4 , 0] U [0, R]. 



4. The GLM representation 

In this section we derive Gelfand-Levitan-Marchenko (GLM) representations for 
the eigenfunctions {3?j}f- These representations will be needed in sectional 



Proposition 4.1. Suppose q(x,t) is a solution of (1.1) with sufficiently smooth 
boundary data go(t) = q{0,t) and gi(t) = q x (0,t), < t < T. The functions $1 
and $2 defined in can be represented as 



$l(t,kj\ _f0\, [' (h(t,s)k + I 3 (t,s)k 3 \ , 

$ 2 (t,k)J - \l) + J_ t \l Q (t, S ) + I 2 (t,s)k 2 J 6 



(4.1) 



where 



I {t,s) 
h(t,s) 



l _ ..ft, ,..,2- --s 1 



j(t,s) + -g (t)l(t,s) + \^-\go(t)\ 2 go(t) + ^gi(t)j n (^ s )^ ( 4 - 2a ) 
l(t,s)- l -g {t)m(t,s), (4.2b) 



I 2 (t, s) = m(t, s) + -go(t)n(t, s) 
h(t,s) = n(t,s), 



(4.2c) 
(4.2d) 



and the functions {n(t, s),m(t, s),l(t, s), j(t, s)} , \s\ < t < T, satisfy the initial con- 
ditions 



m(t,-t)=j{t,-t) = 0, n(t,t)=g {t), l(t, t) = - \g (t)\ 2 g {t) + l -g x (t), 



(4.3) 



and the ODEs 



/ n t (t,s) - n s (t,s) \ 
m t (t,s) +m s (t,s) 
lt(t,s) - l s (t,s) 
\ jt(t,s)+j s (t,s) J 



f2iA 2 (t) ai(t) 

a 2 (t) ax(t) 

a 3 (t) a 4 (i) 2iA 2 (t) a 5 (t) 

\ a 6 (t) a 3 (t) a 4 (i) J 



2fld(*)\ fn(t,s)\ 
m(t, s 
l(t,s) 



(4.4) 
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with the functions {ctj(t)}f given by 
M 2 5o + 2%) , 
|ffo| 4 5o + i§o9i - 2i 9ot) , 

bol 6 + gogot - \go\ 2 go9i + \go\ 2 gogi + *M 2 ) , 

3|5o| 4 5o - 2ig 2 l g 1 + 2i\go\ 2 gi + <%ot) , 
2 7o + 

7|5o| 6 5o - 4ig 2 got - l2i\g \ 2 got 

+ 6i|5o| 2 5oS'i - 12i|#o| 4 5i + 8^0 1 S-i 1 2 - 4ffo5i 
Proof. The matrix-valued function fc) defined by 

<t>(t,k) = fi 2 (0,t,k)e- 2ikHa3 

satisfies 

& + 2ik i a 3 (j) = V 2 <f>, 0(0, fc) = /. 
We seek a GLM representation for (j> of the form 

/t 
[J(t, s) + kL(t, s) + k 2 M(t, s) + k 3 N(t, s)] e- 2ik4sa3 ds, 

(4-7) 

where L, M, N, J are 2 x 2-matrix valued functions of (t, s), \s\ < t < T. Substitution 
of (|4.7[) into (4.6 ) gives 







Ct\ = 


-(1 


2 u 




1 „ 


Ot 2 = 


4 




1 , 


a 3 = 








«4 = 


*(' 


«5 = 


Isol 






«6 = 


32^ 



lit 



(4.5) 
(4.6) 



[J + kL + /c 2 M + k 3 N] e~ + 



J + kL + k l M + k 6 N 

t 



(4.8) 



+ J [J t + kL t + k 2 M t + k 3 N t ] e-ds + 2ik 4 a 3 J [J + kL + k 2 M + k 3 N] e~ds 
= V 2 (e~ + J [J + kL + k 2 M + k 3 N}e~ds^ , 



where we use the following short-hand notations: 



-2ik 4 ta 3 



e := e 



and 



F := F(t,t), F := F(t, -t) whenever F = N, M, L, or J; 

all other occurences of the functions N, M, L, J and their derivatives are understood 
to be evaluated at (t, s); if the exponential e~ is part of an integrand, it is understood 
to be evaluated at s instead of at t. 



We substitute the expression (2.5) for V 2 into equation (4.8) and apply the inte- 
gration by parts identity 



-2ik 4 



Fe~ds 



Fe~ - Fe 



F^e ds 



0- 3 



to all terms in the resulting equation which are of order 0(k J ) with j > 4. This 
yields an equation only containing terms of order 0{k 3 ) with < j < 3. Equating 
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the coefficients of the terms of order k 3 which do not contain integrals, we find 

(N - a 3 Na 3 )e' + (N + a 3 Na 3 )e = 2Qe~ . 

Equating the coefficients of the terms of order 0(k 3 ) which do contain integrals, we 
find 

/N t e~ds + a 3 I N s e~dsa 3 = 2Q / Je~ ds — iQ 2 cr 3 / Le~ ds 
-t J-t J-t J-t 

+ (Q 3 - iQicrs) [ Me~ds - iA 2 a 3 f Ne~ds 
J-t J-t 

We deduce that the terms of order 0(k 3 ) vanish provided that 
N - a 3 Na 3 = 2Q, N + a 3 Na 3 = 0, 

N t + a 3 N s a 3 = 2QJ - iQ 2 a 3 L + (Q 3 - i&as) M - iA 2 a 3 N. 

Similarly, the terms of 0(k 2 ) vanish provided that 

M - <r 3 Mo 3 = -iQ 2 a 3 + iQNa 3 , M + a 3 Ma 3 = -iQNa 3 , 
M t + a 3 M s a 3 = -iQ 2 a 3 J + (Q 3 - iQ x a 3 ) L - iQN s a 3 - iA 2 a 3 M; 

the terms of 0{k) vanish provided that 

L - a 3 La 3 = Q 3 - iQi<j 3 + iQMa 3 + -Q 2 a 3 Na 3 , 

1 . 

L + a 3 La 3 = -iQMa 3 - -Q z a 3 Na 3 , 

Lt + <t 3 L s cj 3 = (Q 3 - iQia 3 ) J - iQM s a 3 - ^Q 2 a 3 N s a 3 - iA 2 a 3 L; 

and the terms of 0(1) vanish provided that 

1 i 

J - a 3 Ja 3 = iQLa 3 + -Q 2 a 3 Ma 3 + -(Q 3 - iQi<r 3 )Na 3 - iA 2 a 3 , 

1 ~ i 

J + a 3 Ja 3 = -iQLa 3 - -Q 2 a 3 Ma 3 - -(Q 3 - iQia 3 )Na 3 , 

1 i 

J t + a 3 J s a 3 = -iQL s a 3 - -Q 2 a 3 M s a 3 - - (Q 3 - iQi<r 3 ) N s a 3 - iA 2 a 3 J. 
These equations suggest that we introduce L, M, J by 
M = M + l -Qu 3 N, 

L = L+ l -Qa 3 M, (4.9) 
J = J+ l -Qa 3 L+ fjQVa + iQi^ N. 



2 V 8 4 ' 

Straightforward, but tedious, computations show that 

F + <j 3 F(t 3 = whenever F = N, M, L, or J, 
N(t,t)-a 3 N{t,t)a 3 = 2Q, 

M(t, t) - a 3 M(t, t)a 3 = 0, (4.10) 
L(t, t) - a 3 L{t, t)a 3 = Q 3 - iQ x a 3 , 
J{t,t)-a 3 J(t,t)a 3 = 0, 
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and 

N t + a 3 N s a 3 

M t + a 3 M s a 3 

L t + a 3 L s a 3 



Q 3 - iQi^s M + 2QJ 



--Q 5 --A 2 Q--Q t a 3 + -Q 2 Q 1 a 3 jN+^-Q 3 -iQ 1 a 3 jL 
l -Q 6 a 3 - % -A 2 Q 2 a 3 + l -QQ t - l -Q 2 Q x Q + \q 3 Qi - ^QiQ 3 



+ - A Q\<J 3 )N+ ( ~-Qtcr 3 - -A 2 Q\M + (Q 3 - iQia 3 ) J, (4.11 



Jt + cr 3 J s a 3 



The symmetries 



- Q(T:i ( L, + a :i L*a 3 ) + ( -^Q 3 <7 3 - ^Ql ) ( A'/ + <T:!.A>:i ) 



2 QtO" 3 L - ^(Q 3 )to- 3 iV - -Qi t JV - iA 2 a 3 J. 



cf)(t, k) = a 3 <p(t, -k)a 3 , 4>{t, k) = <7i0(i, fc)oi, 

suggest that we seek L, M, N, J such that J and M (resp. L and N) are diagonal 
(resp. off-diagonal) matrices, and such that 

F = a\Fa\ whenever F = N, M, L, or J. 

In fact, introducing complex-valued functions {n(t, s), m(t, s),l(t, s),j(t, s)} by 



N 



L 





ne J (°' ) 



ne J (°' ' 



r(o,t) 



M 



m 
m / ' 



/e 2lJ (o, 0) A 

2i f (0,t) A 

le J (°.°) 



J 



J 



(4.12) 



a long computation shows that the two sets of equations (4.10) and (4.11 ) are satisfied 
provided that {n,m,l, j} satisfy the initial conditions (4.3) and the ODEs (4.4), 
respectively. 



Equations (|4.7|), (|4.9|), and (|4.12|) imply that 

rt ( h(t, s)k + I 3 (t, s)k 3 \ 2ik i s 



e 2iJ (°$ A <t> 12 (t,k) 

4>22{t,k) 





,2ikH I + 



I (t,s) + I 2 (t,s)k 2 



2 e 2lk s ds, 



where {Ij}^ are given by (4.2). The relations 

12 = d> ie -2^o tA 2 e 2ifc 4 t ) 

then imply that {^j}? are given by (4.1). 



■>22 



$2e 



2ik 4 t 



□ 



5. The solution of the global relation: physical domain 

In this section, we consider the solution of the global relation in the physical do- 
main. Theorems 5.1 and 5.3 below lead to expressions for the Dirichlet-to-Neumann 



map in terms of the solution of a system of nonlinear integral equations. Theorem 5.1 



applies to arbitrary initial data. Theorem |5.3| provides a more convenient expression 
for the Dirichlet-to-Neumann map in the special, but important, case of vanishing 
initial data. 
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5.1. The general case. For a function F(t, s), — t < s < t, we let A denote the Abel 
transform in the second variable defined by 

(AF)(t,s)= / F(t,s>) 



j-t Vs - s' 

The inverse A -1 of this transform is given by 



(A-'F^s) 



1 d f s F(t, s')ds' 1 



7T ds 



t Vs — s' 



7T 



l£(t,s')d S ' , F(t,-t) 



+ 



s — s' 



Vt + s 



Theorem 5.1. Let qo(x), go(t), a(k), andb(k) be as in theorem 3.1 and let {kj, —kj}^ C 
D\ U Z?2 denote the possible simple zeros of a(k). Define functions fo(t) and f2(t), 
t > 0, by 



7' 



j = 0, 2, t > 0, 



(5.1) 



where 7' is f/ie contour defined in section^ Let q(x,t), x > 0, < i < T, be a 
solution of the DNLS equation (1.1) which satisfies 

q(x,0) =q (x), x>0; q(0,t) = g {t), < t < T. 



Let {n,m,l,j} be the functions in the GLM representation (4-1)- 
Then the Neumann boundary data gi(t) := q x (0,t) are given by 



= i\go{t)\ 2 g {t) - 4ci(.A- 1 n)(M) + go(t)m(t,t) + 



(o,t) 

4 e "'J(o,o) ' 



7T 



2*/ 2 (2t) 



+ 2i 



2i 



- ( m(i,i) - glfloC*)! ) /o( 2 *) + Ygflo(t)n(t, -t)?o(0) 



f2(t + s)ds 
(5.2) 



+ 



m s(i, s ) - ^go{~t) n s{t, s) ) f (t + s)ds 



where c\ E C is a constant defined by 

ci = 



le~ 2u4 dl 



ire 



-J7T/4 



2V2 



(5.3) 



Proof. Substituting the GLM representation (4.1) into the global relation (2.13a), we 
find the following equation which is valid for k E Dx U F>2'- 



-t 



{h(t,s')k + I 3 (t,s')k 3 ) e **V-*W 



a(fc) 



1 + y" (l (t, s') + J 2 (t, S ')& 2 ) e- 2ife4 ( s '-*)d S ' 



We multiply this equation by k 2 e 2ik4 ^~ s \ —t<s<t, and integrate the resulting 
equation with respect to dk along 7'. This yields 



(h(t, s')k + I 3 (t, s')k 3 ) e 2ik ^ s '-^ds'dk + eHo/o) A 



M) 



a{k) 



(5.4) 
dfc = 0, 
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where the right-hand side is zero because of Jordan's lemma. Analyticity implies that 



the contour 7' can be replaced by 7 in the first integral in (5.4). Using the Fourier 
inversion identitjj^] 



J k 3 J* K(t,s')e 2ik ^ s '- s Us' 



dk = \i K ^ s )i -t<s<t, 

ifiT(M), 0<s = t, 



we see that the term in (5.4) involving I\ equals 

ft 



[ k 3 [ I 1 (t,s')e 2ik ^ s '- s Us'dk = ^I 1 (t,s). 
Jy J-t 4 



(5.5) 



(5.6) 



In order to simplify the term in (5.4) involving ^3, we integrate by parts as follows: 



k 5 f h(t,s')e 2iki{s '^ds'dk 
7 J-t 

1 J k /3(^,^)e 2^fc4 ( ^ - s )-/3(^,-^)e- 2^fc4(^+s) -y'%3s'^^ , )e 2^fc4(s '- s) ^ , 



dk. 



We can now interchange the order of integration in the integral on the right-hand 
side involving /3 S /. Using the identity 



(5.7) 



where c\ is given by ( |5.3[ ), we find that the term in (5.4) involving ^3 equals 

a da! 



2i 



h{t,-t) r , 

yt + s J-t v s - s 



ciir 



{A- l h){t,s). (5. 



In view of the definition (|5.1|) of /o and f 2 , the term in (5.4) involving Jo equals 

I (t,s')f 2 (s + s')ds'. (5.9) 



In order to simplify the term in (5.4) involving I 2 , we integrate by parts as follows: 



fc 4 ^ f I 2 (t,s')e- 2ik4(s ' +s) ds'dk 
r a[k) J_ t 



2i Jy a(k) 



h(t,t)e 



-2ik 4 (t+s) 



h(t,-t)e 



I 2 At,s')e- 2ik ^ s ' +s Us' 



-2ik 4 (-t+s) 



dk. 



Expressing the integrals with respect to dk on the right-hand side in terms of /o, we 
infer that the term in (5.4) involving I 2 equals 



21 r ( ° !t) a 

2i 



h(t,t)f (t + s) - I 2 (t, -t)f (-t + s) - / l2s'(t,s')f (s' + s)ds' 



(5.10) 



2 This identity can be less rigorously expressed as J fe 3 e 2lfc ^ s ^dk = 
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Using equations (5.6) and (5.8)-(5.10), we can write (5.4) as 



1 

2d 



h(t, t)/o(* + s)- I 2 (t, -t)f (-t + a) 



hs'(t,s')f {s' + s)ds' 



-t 



0, 



Recalling the definitions (4.2) of {/j}o and the initial conditions (4.3), evaluation of 
this equation as s — > t yields 



l(\\go(t)\ 2 go(t) + % -gi{t) - l - go (t)m(t,t) 



+ 



1 



j(t,s) - -g (t)l(t,s) + -\g (t)\ 2 g Q (t) + -l(t,t) n(t,s) 



2i 



f 2 (t + s)ds 



1 

2% 



m(t,t) - -\g (t)\ 2 / (2t) + -g (t)n(t,-t)fo(0) 



(5.11) 



m s (t, s) - ~go(t)n a (t, s) f Q (t + s)ds 



Since 



we find 



ft(fc) ^o(O) 
a(jfej 



2fc 



O 



oo, /ceDiUD 2 , 



/o(0) 



7' 



a{k) 



dk 



7' 



»go(0) 
2k 



dk 



7rgo(0) 



Solving equation (5.11) for fl'i(t) and using (5.12), we find (5.2). 



(5.12) 

□ 



Substitution of the expression \h.2\ for gi{t) into (4.4) yields a nonlinear system of 
integro-differential equations involving {n,m,l, j} and A2. Employing the identities 



(d t -d s )F(t,s) = f(t,s) 



(d t + d s )F(t,s) = f(t,s) 



F(t,s)=F 



F(t,s) = F 



t + s t + s 



t — s t — s 



+ I f (r,t + s - r) dr, 

Jt±s 

(5.13) 

+ / /(r,r -t + s)dr, 



t-s 
2 



and recalling the initial conditions in (4.3), we can convert this into a system of 
nonlinear integral equations. Supplementing this system with the equation for A2 

the exp 

A 2 (0,t) 



obtained by substituting the expression (5.2) for gi(t) into the equation 

3 



Ifol 



:(gi9o ~ 9o9l), 



4 IJU| 2 

we find a closed system of nonlinear integral equations for the complex-valued func- 
tions {n(t, s),m(t, s), l(t, s),j(t, s)}, \s\ <t<T and the real- valued function Aa(0, t), 
< t < T. This system is formulated explicitly in terms of the initial data qo(x) and 
the Dirichlet boundary data go(t)- Once the solution of this system has been found j^] 
the Neumann data gi(t) can be determined from (5.2). This provides an effective 
way of computing the Dirichlet-to-Neumann map. 



The problem of rigorously showing existence and uniqueness for this system remains open. 
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IV 



Remark 5.2. The expression (5.2) for g± in terms of {n,m,l, j} can also be derived 

states 

(5.14) 



directly from the expression (3.3) for g± given in theorem 3.1. Indeed, (3.3) states 
that g\ is given by 

*l 12 
2 



91 



\g \~g + 4c {3) + 2ig ^> 



(2) 
2 i 



where is given by (3.4) and 



(2) 



2 



k[$ 2 (t,k) - <f> 2 (t,ik)]dk. 



(5.15) 



Substituting the GLM representation (4.1) for $2 into (5.15) and using the identity 
(5.5), we conclude that 



TT 



h{t,s)e 2ik4 ' {s ~ t) dsdk = g (m(t,t) + ^gon(t,t)) . 



On the other hand, we can write as the sum of two terms, 



-(3) (3) 
•T ~r L 2 ' 



,(3) 



k 2 $x(t,k) + i$i(t,ik) + 



k 



dk, 



where 



(5.16) 



(3) 

and Cn in volves the quotient b(k)/a(k). Substituting the GLM representation for $1 
into (5.16) and integrating the term in the resulting expression involving A; 5 by parts, 

-AikH 



we find 



.(3) 



kn(t, t) — kn(t, —t)e 



j kn s (t,s)e 2lk4 ^d& 



kgc 



dk. 



Using the identity (5.7) and the initial condition n(t,t) = go(t), this gives 



.(3) 



-ci(A~ l n)(t,t). 



(5.17) 



Finally, simplifying c)> as in the proof of theorem 5.1 we deduce that c^' is given by 



.(3) 



the terms on the right-hand side of equation (5.2) which involve the functions /o and 
f%. In summary, after collecting the various contributions, we conclude that (5.2) can 
be derived from (3.3). 

5.2. The case of vanishing initial data. With respect to the integral equations 
derived in the spectral domain (see theorem 3.1), the system of integral equations 
derived in section 5.1 has the numerical advantage of being defined on the bounded 
domain |s| < t < T. Nevertheless, the system of section 5.1 involves five highly 
nonlinear and rather intricate equations. The next theorem shows that this system 
can be simplified considerably in the case of vanishing initial data qo(x) = 0. 

Theorem 5.3. Suppose that qo(x) = 0, x > 0. Let go(t) and q(x,t) be as in theorem 



5.1. Then the Neumann boundary data gi(t) := q x (0,t) are given by 
9i(t) = i\9o{t)\ 2 g {t) - 4ci(^ 1 n)(t,t) + g (t)m(t,t), 



(5.18) 



following system of two nonlinear integral equations: 



where c\ is given by (5.3) and the functions {n(t, s),m(t, s)} , \s\ < t < T, satisfy the 

(5.19a) 



nit, s) 



m(t, s) 



t + s 

go l — — I + 



t+s 
2 



[Pi(r)n + /3 2 (r)m + P^A^m)} dr, 



t -s 
2 



[Pi(r)m + P 5 {r)n + ^{A^n)] dr, 



(5.19b) 
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5i 1 

j |50 1 4 + ^^oSi - 5051, 



= \go\ go + igi, 
'2ci' 



ft 
/3 4 



1 

4 
1 



Iffol 4 + 2#o5i) , 



(5.20) 



ft = j [\ga\ 4 go + ^o5i - zigot) , 



ft 



7T 

8c^ 



|5o| 2 5o - 2%) • 



Unless otherwise specified, the functions in (5.19a\ ) are evaluated at (r,t + s — r) and 
£/ie functions in (5.19b) are evaluated at (r, r — i + s). 



Proof. The expression (5.18) for gi(i) is obtained by setting fo(k) = f2(k) = <7o(0) = 
in (5.2). We will use the global relation to eliminate the functions I and j from the 
equations (|4.4[) for {n,m,Z, j}; this will lead to (5.19). 



Substituting the GLM representations (4.1) into the global relation (2.13) and 
using that b(k) = 0, we find 



y* (fell + £; 3 /3)e 2ifc4 ( s '-'W = c(Jfe, i), 
y* (I + I 2 k 2 )e 2ik ^ s '~^ds' = d(k, t) - 1, 



(5.21a) 



(5.21b) 



where c and d — 1 are analytic functions of £ Z?i U D2 and of order 0(l/k) as 
— )• 00. 

We first analyze (5.21a). We multiply (5.21a) by e 2lk (*~ s ) with —t<s<t and 



integrate the resulting equation with respect to dk along 7. This yields 

{kh + k 3 I 3 )(t, s')e 2iki{s '- s) ds'dk = 0, 



(5.22) 



■yJ-t 



where the contribution of the right-hand side vanishes by Jordan's lemma. Using the 
identities (5.5) and (5.7), equation (5.22) becomes 



7T 



Cl (Ah){t,s) + -I 3 {t,s) = 0, 



-t < s <t. 



(5.23) 



Applying t he in verse Abel transform to (|5.23|) and using the expressions for I± and 
I3 given in (4.2), we fine 4 ' 



KM) - 7,go{t)m(t,s) 



+ ^{A' 1 n)(t,s) =0, -t<s<t. 



(5.24) 



Similarly, we multiply (5.21b) by ke 2tk (* s ) with —t<s<t and integrate the 
resulting equation with respect to dk along 7. Steps analogous to those that led to 

(5.25) 



(5.24) now yield 
ci 



1 



j(t,s) + -g {t)l(t,s) + -\g (t)\ 2 g (t) + 7 gi(t) n(t,s) 



+ jiA^mXt, s) + ySoCtX^nXt, s) = 0, 



-t < s <t. 



4 In view of (4.3|, evaluation of (5.24| as s — > t yields an alternative derivation of (5.181. 
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We solve the equations (5.24) and (5.25) for j and I and plug the result into the first 



two equations of (4.4). This yields 



n t - n s 
nrit + m< 



Pi 

ps + PsA~ l 



+ PsA' 1 
Pa 



where {Pj}\ are given by (5.20). Using the identities (5.13) and the initial conditions 
(Ob, we find d57l9b. □ 



Substitution of the expression (5.18) for gi(t) into (5.19) yields a system of nonlin- 
ear integral equations for {n(t,s),m(t,s)}, \s\ < t < T, formulated entirely in terms 
of the known boundary values go{t). This system only involves two equations for two 
unknown functions. This is in contrast to the system derived in section 5.1 in the 
case of general initial data, which involves five equations. 
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